Chaotic Motion Around Prolate Deformed Bodies 
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The motion of particles in the field of forces associated to an axially symmetric attraction center 
modeled by a monopolar term plus a prolate quadrupole deformation are studied using Poincare 
surface of sections and Lyapunov characteristic numbers. We find chaotic motion for certain values 
of the parameters, and that the instability of the orbits increases when the quadrupole parameter 
increases. A general relativistic analogue is briefly discussed. 
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Attraction forces represented by a monopolar plus a 
prolate quadrupolar distribution of masses (charges) are 
a good approximation for elongated massive (charged) 
bodies. Examples range from astrophysics to nuclear 
physics. There are many observed galaxy clusters with a 
cigar like shape H . Also, the nuclear charge of light gold 
atoms has been reported as having a large prolate defor- 
mation 0. Most of the Dwarf Galaxies in the Virgo 
Cluster may obey the "prolate hypothesis", i.e., they 
probably have a prolate spheroidal shape ||. Asteroids 
also have a prolate shape, but usually they are not ax- 
isymmetric. Merrit Q found, from detailed modeling of 
triaxial galaxies, that most of the galaxies must be nearly 
axisymmetric, either prolate or oblate. 

Classical, as well as, quantum chaos have been studied 
in a variety of axially symmetric fields of forces. In partic- 
ular, attraction centers described by potentials that are 
the sum of two terms: a monopolar term and a quadrupo- 
lar deformation. Furthermore this center is "perturbed" 
by an external distribution of masses (charges) repre- 
sented by its external multipolar moments, i.e., 

V= -a/R-qP 2 {costf)/R 2 + V>, (1) 
V P = Qii?Pi(cosz?) + Q 2 R 2 P 2 (cos &) + ... . (2) 

Sometimes the monopolar term is changed by the poten- 
tial of a spring In general, in all these cases the terms 
that originate the chaos are the external multipolar mo- 
ments. 

We shall consider the simplest, albeit, important case 
of a particle moving in the field of a monopolc plus a 
quadrupole deformation. This deformation is usually 
considered to be the major deviation from spherical sym- 
metry. In cylindrical coordinates, (r, tp, z), the field takes 
the generic form, 

tji ^ g g(2z 2 -r 2 ) 

U ( r > z ) = / 9 , 9 ~ of 2 i 2A5/2 ' ( 3 ) 

Vr 2 + z 2 2(H + z 1 )*! 2 - 

where a is a constant that may be associated with the 
central body mass (charge). It is instructive to have a 
special model in mind, consider two equal masses located 
on the z-axis symmetrically, at z = —a and z — +a. The 
gravitational potential of the above mass configuration 



up to the order a 3 is (||) with q — lac? . We shall use 
a = 1 without loss of generality. Note that we are not 
considering external multipolar moments (Vp = 0), i.e., 
only deformed cores will be studied. 

We can distinguish two cases depending on the sign of 
q. The oblate deformation case, q < 0. This is the com- 
mon case for bodies deformed by rotation and has been 
analyzed in astronomy for more than two hundred years. 
The integrability of the Newton equations for a particle 
moving in the gravitational field of an axially symmetric 
oblate body is an unsolved problem. It is known as the 
classical problem of the existence of the third isolating in- 
tegral of motion . There are numerical evidences that 
orbits of particles moving around a monopole plus an 
oblate quadrupole are not chaotic. We study the prolate 
deformation case, q > 0; that is the one that we shall 
discuss in this communication. In this case we have a 
monopolar field (the usual Kepler problem) "perturbed" 
by a quadrupolar term, in other words, we have a typical 
situation wherein the KAM (Kolmogorov-Arnold-Moser) 
theory applies 0. 

First we study the contours of the effective potential 
U e ff = U + h 2 /(2r 2 ), where h z — r 2 ip is the axial specific 
angular momentum that due to the axial symmetry is 
conserved. We also have the conservation of the total 
specific energy, E = (f 2 +z 2 )/2+U e ff. Thus we have that 
the motion is completely determined by the functions 
r = r{t) and z — z(t). Then, we have a four dimensional 
phase space. But, due to energy conservation the motion 
actually takes place in a three dimensional space. An 
adequate tool to investigate the trajectories in this phase 
space is the Poincare surface of section method. Now 
let us comeback to the effective potential contours. In 
Fig. 1 we plot the level contours of U e fj for L z = 0.83 
and E = 0.464 and different values of the quadrupole 
moment parameter: a) q = 0.3, b) q = 0.5, c) q = 0.85, 
and d) q=0.95. Thus for these values of the parameters 
the motion of the particle is confined to toroidal regions 
that do not contain the symmetry axis. Note that for the 
last case we have two non connected regions. 

The particles move in the reduced phase space (p r = 
r,r,z). Note that p z — z is determined by the energy 
conservation. In Fig. 2, for the case a), we present the 
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intersection points of some particle trajectories with the 
plane z = 0. The picture is the one for regular orbits. 
The case b) is analyzed in Fig. 3, using the same surface 
section. We find regions of non destroyed tori together 
with chaotic regions in concordance with the KAM the- 
ory. In Fig. 4 we show again the case b) but, now with 
a different section, z = 0.4. We see that the integrable 
and chaotic regions are deformed depending on the cho- 
sen section. We also studied the case c) that is quite 
similar to the former, so we shall not present it here. We 
find that increasing the quadrupole moment the size of 
the chaotic regions also increases. And finally, in Fig. 5 
we study orbits in one of the non connected regions of 
the case d). In this last case the surface section is taken 
as z — 0.4, again we find large regions of chaotic behav- 
ior and some non destroyed tori. In summary, we find 
chaotic behavior of orbits for several values of prolate 
quadrupole moment. 

To quantify the degree of instability of the orbits 
we shall study their associated Lyapunov characteristic 
numbers (LCN) that are defined as the double limit 



LCN = lim 

<5o -> 



log(J/(Jo) 



(4) 



where 5o and S are the deviation of two nearby orbits at 
times and t respectively. We get the largest LCN by 
using the technique suggested by Benettin et al. 

We fix the value of the constants of motion as L z = 
0.83 and E = 0.464, and choose the same values of 
quadrupole parameters used to plot the Poincare sec- 
tions. For the value q — 0.3 it was chosen the refer- 
ence orbit with initial conditions: z = 0, p r — 0, and 
r = 0.85, and S ~ 10~ 9 , we found LCN < 10~ 4 that 
characterizes a stable system. With q = 0.5 and initial 
conditions: z = 0, p r = 0.2, Finally, for q — 0.95, and 
z = 0.4 , Pr = 0.05, and r = 0.95, we obtain LCN ~ 0.09 
(±0.015). We see that the degree of instability increases 
when the quadrupole parameter increases for fixed con- 
stants of motion. 

As we said before, there are numerical evidences that 
orbits of particles moving around a monopole plus an 
oblate quadrupole are not chaotic. The difference be- 
tween the oblate and the prolate case can be under- 
stood analyzing the critical points of the effective po- 
tential U e ff- In particular, the existence of the sad- 
dle points that is one of the main ingredients to have 
instable motion. We find that the critical point, r — 
a/ (3q + 2L|) /2a, z — 0, is a saddle point if the pa- 
rameters obey the two conditions, L 2 < 3q and 3L 2 > 
y2o/3. therefore when q < 0, the oblate case, no real 
L z can obey the first of these condition. 

The Newtonian motion of a particle moving in the po- 
tential (3) has a general relativistic analogue. The po- 
tential is replaced by a metric solution to the vacuum 
Einstein equation and the particle motion equation by 



the geodesic equation. The instability of geodesies in 
metrics associated to a black hole surrounded by a shell 
of matter was studied in some detail in Q . 

A solution to the Einstein equations that has as a 
Newtonian limit a potential like (||) is the Erez-Rosen- 
Quevedo (ERQ) solution We did not found chaos 

in the oblate case, the prolate case is also chaotic. The 
confinement region for the relativistic motion constants, 
E = 0.937, and L z — 3.322, and the quadrupole param- 
eter q — 5.02 is presented in Fig. 6. The coordinates used 
in this case are spheroidal coordinates, that are the ones 
appropriate for the ERQ solution, they are related to the 
usual cylindrical coordinates by u = + B_)/(2m), 
and v = (R + -R_)/(2m), with R ± = (V 2 + (z±m) 2 ) 1 / 2 . 
We have two regions of confinements that we have la- 
beled II and III. In Fig. 7 we present a Poincare section 
for particles moving in the region II, the section is taken 
as v = 0, i.e., z = 0; u,p u are canonical conjugate vari- 
ables. We see a phase space with a large region of chaotic 
motion. We shall present a complete study of geodesic 
in ERQ spacetimes elsewhere. 

Some dense cores in dark clouds have been found to 
have prolate spheroidal shape [Q . Then a prolate geom- 
etry has to be considered as initial condition in the star 
formation process. We think that the strong instability 
presented here may play a crucial effect in the formation 
of structures in starts O] . 

We want to finish this short communication by re- 
minding that in nonlinear systems of equations chaos is 
the rule rather than the exception. Thus simple systems 
with a minimum of structure play an important role in 
the physical, as well as, mathematical understanding of 
chaos. A good example is the paradigmatic Henon-Hciles 
system wherein the "simple" addition of a x 2 y term in the 
potential of two uncoupled oscillators (integrable motion) 
has dramatic consequences that are the physical manifes- 
tation of the creation of a saddle point together with a 
perturbation. In the case presented here, we have a sim- 
ilar situation, the prolate-quadrupole potential also adds 
a saddle and a perturbation. 
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FIG. 1. We plot the level contours of U eff for L z = 0.83 
and E = 0.464 and different values of the quadrupole moment 
parameter: a) q = 0.3, b) q = 0.5, c) q = 0.85, and d) q=0.95. 
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FIG. 2. Surface of section for L z = 0.83 and E = 0.464 
and q = 0.3. The section corresponds to the plane z — 0. For 
these values of the parameters we have the section of regular 
motion. 



0.4 




' 0.65 0.85 1 .05 1 .25 1 .45 1 .65 

r 

FIG. 3. Surface of section for L z = 0.83 and E = 0.464 
and q = 0.5. The section corresponds to the plane z = 0. For 
these values of the parameters we have the typical section 
indicating chaotic motion. 
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FIG. 4. Surface of section for the same values of the pa- 
rameters that in the precedent figure, but a different section, 
z — 0.4. We see a different cut of the regular and chaotic 
regions. 
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FIG. 5. Surface of section for L z = 0.83 and E = 0.464 
and q = 0.95. The section corresponds to the plane z = 0.4. 
Again we have irregular motion. 
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FIG. 6. Level contour for the general relativistic 
quadrupole + monopole system (ERQ solution). The rela- 
tivistic constants are L z = 3.32 and E = 0.937 and q = 5.02. 
The labels u and v denote spheroidal coordinates. 



FIG. 7. Surface of section for the region II shown in pre- 
vious figure. We have a large region of chaotic motion. The 
section corresponds to the plane v = i.e., 2 = 0. 
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